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We study the consequences of random mass, random scalar potential and random vector potential 

on the line of fixed points between integer/fractional quantum Hall states and an insulator. This line 

of fixed points was first identified in a clean Dirac fermion system with both Chern-Simon coupling 

and Coulomb interaction in Phys. Rev. Lett. 80, 5409 (1998). By performing a Renormalization 

Group analysis in 1/A^ ( A'' is the No. of species of Dirac fermions) and the variances of three 

disorders Am, Av, Aa, we find that Aa/ is irrelevant along this line, both Aa and Ay are marginal. 

With the presence of all the three disorders, the pure fixed line is unstable. Setting Chern-Simon 

interaction to zero, we find one non-trivial lines of fixed points in {/S.a,w) plane with dynamic 

^v' ' exponent z = 1 and continuously changing v, it is stable against small (Am, Av) in a small range 

Q>^ I of the line 1 < w < 1.31, therefore it may be relevant to integer quantum Hall transition. Setting 

Aa/ = 0, we find a fixed plane with z = 1, the part of this plane with i^ > 1 is stable against small 

-^ Aa/ , therefore it may be relevant to fractional quantum Hall transition. Although we do not find a 

' <j -J , generic fixed point with all the couplings non-vanishing, we prove that the theory is renormalizable to 

fT , ' the order (1/7V)^, (1/A'^)A, A^ and explore the interesting processes which describe the interferences 

between Chern-Simon interaction. Coulomb interaction and the three kinds of disorders. 
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^^ [ The zero temperature quantum phase transitions between the different quantum Hall and insulating states of a two- 

O^ I dimensional electijee gas in a strong magnetic field are among the most intensively studied quantum critical points, 
i^ ■ both theoreticallyoa and experimentally^. Earlier theoretical investigations focussed on the transitions between the 
Integer Quantum Hall plateaus and described them in terms of non-interacting electrons moving in a random external 
potentialQ. 

""^ , Ludwig et. al introduced and analyzed a Dirac fermion model with random mass, random scalar potential and ran- 

i-rt ' dom gauge potential. They found that the random mass is marginally irrelevant, random scalar potential marginally 
P3 I relevant. However, they found that random vector potential is exactly marginal, therefore there is a line of fixed points 
O ■ characterized by the strength of the random gauge potential- the zero energy wavefunction shows multifractal behav- 
^ ' iors with the exponents continuously changing along the lingj. If two of the random potentials are non-zero, the third 
^ one will be generated and the system flows to strong coupling regime. They further argued that the system should flow 
to thej-generic fixed point of integer Quantum Hall transitions which was correctly described by Chalker-Coddington 
modeiotl. 
j^ ' The|PCoperties of the fixed line of Dirac fermions in the presence of random magnetic fields-sYere further investigated 

in Ref.&Q. Later, the U{1) gauge potential was extended to non-Abelian gauge potentia]li3~E3. 

It has also been argued that the transitions between fcactional quantum Hall states could be mapped onto models 
essentially equivalent to those betwecp^thc integer statest^. The latter point of view was however questioned by Wen 
and WuEJ and Chen, Fisher and WullJ: they focussed on the simpler case of systems in the presence of a periodic 
rather than a random potential, and examined a model of anyons, with a statistical angle 9 and short-range repulsive 
interactions, which displayed a second order quantum phase transition between a quantized Hall state and a Mott 
insulator as the strength of the periodic potential was varied. This transition was characterized by a line of critical 
points with continuously varying exponents, parameterized by the value of 9. For the case 9 = 0, when the anyons 
were fermions, the transition was out of a integer quantum Hall state; its exponents and other universal properties 
were different from the cases < 6* < 27r for which the anyons acquired fractional statistics and the transition was out 
a fractional quantum Hall state. (For 9 ^ 2tt the anyons became bosons and the Hall state reduced to a superfiuid.) 
In all of the above theoretical works, the long-range Coulomb interactions between charge carriers have beea 
effectively ignored However, a few recent works have taken steps to remedy this serious shortcoming. Yang et alx3 
studied] the integer quantum Hall transition under a Hatree-Fock treatment of the Coulomb interaction. Lee and 
WangllZI showed that the renormalization group eigenvalue of the Coulomb interaction was zero at the Hartree- 
Fock criticaLpoint; higher order calculations are therefore necessary to understand the physics. Pfannkuche and 
MacDonaldli3 numerically studied electrons with Coulomb interactions in a periodic potential between a fractional 



Hall state and an insulator, but were limited to rather small system sizes. Interesting scaling interpretations of 
Coulomb interaction-induced dephasing were discussed in ReftJ. 

Most recently, neglecting the disorders. Ye and Sachdev provided-a-Uiorough analysis of the consequences of Coulomb 
interactions on the anyons in a periodic potential model of Refgy'tJ. They showed that the Coulomb interaction is 
marginally irrelevant for the integer case (0 — 0), and remains so for the fractional case for small values of 9; this 
marginally irrelevant interaction will lead to logarithmic corrections to naive scaling functions for the vicinity of the 
transition. For larger 0, they established, in a certain 1/A'^ expansion, the existence of a novel line of fixed points 
at which the Coulomb interactions acquire a non-zero fixed point value determined by the value of 9. There are no 
logarithmic corrections at these fixed points, and naive scaling holds. Theyiound a dynamic critical exponent z — 1 
at all points on the fixed line, providing a concrete realization of the scenariooB, not previously established explicitly, 
that energies must scale as inverse distances for the l/r Coulomb interaction. They also found the correlation length 
exponent v > 2/d ( where d is the spatial dimensionality ) along this fixed line, which implies that the fixed line is 
stable against weak random mass disorder. 

In this paper, we study in detail the consequences of random mass, random scalar potential and random vector 
potential at the line of fixed points. In the rest of this section, we introduce the notations and the model. In the next 
section, we give very general Renormalization Group (RGO formulation of the model and establish some exact results 
( for example, Eqs.|5| and g ). In Sec. Ill, we perform one-loop expansion and discuss the implications of R.G. flow 
equations. In Sec. IV, we flrst discuss the R.G. equation at A^ = oo limit, then calculate 1/iV corrections and discuss 
the solutions in different cases. Finally, we reach conclusions in Sec. V. In the appendix, we show the equivalence of 
the two forms of the random vector gauge potential. |_. 

We begin our analysis by writing down the model of Re£3 extended to include random mass M(x), random scalar 
potential V{x) and random vector potential Ai{xp. 



S = / d'^xdr 



V jV y N 
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The V'm are m — \ ■ ■ ■ N species of charge q/yN 2+1 dimensional Dirac fermions which interact with a U{1) gauge 
field (ao, Qi) {i = 1, 2); we are interested in the case A^ = 1 but will find the large N expansion to be a useful tool. 
The 7o,7i are the Dirac 7 matrices, Xi (r) are spatial (temporal) co-ordinates with d^ = dr, di = dxi, and fc, uj 
(k = \k\) are the Fourier transformed wavevector and frequency variables. To aid the subsequent renormalization 
group analysis, we are working in d = 2 + e spatial dimensions and /i is a renarmalization scale. The parameter a 
is introduced to allow for anisotropic renormalization between space and timg^. We have used the Coulomb gauge 
which allows us to explicitly represent a, in terms of the transverse spatial component with a^ = icijkjat/k. The a^at 
term in S is the Chern Simons coupling: it turns the Dirac particles into anyons with a statistical angle 6 /N with 
= qg; notice that the angle is of order 1/iV and so the expected periodicity of the physics under 6/N -^ 9/N + An 
will not be visible ia the 1/N expansion. The atat term is the Coulomb interaction, and it has been written in terms 
of at following Re£3. ■— ■ 

In the absence of the Coulomb interaction and the random potential terms, it was shown in Reill3 that S represents 
the critical theory of a system of anyons in a periodic potential undergoing a transition from an insulator with 
conductivities Gxx = o'xy — into a fractional quantum Hall state with a^x =^ and a^y = {q^ /h)/{l — 9/27:) to 
leading order in 1/iV. Both these states have energy gaps. It was shown in Ref.EHl that the Coulomb interaction and 
disorder do not modify the values of cjij in either phase. I— . 

The relationship of the continuum model S to the more realistic model of electrons studied in Reio remains 
somewhat unclear, although it is plausible that S is the critical theory of the latter. We may also view S as the 
simplest theory consistent with the following requirements, and therefore worthy of further study: (i) the two phases 
on either side of the critical point have the correct values of Uij , and the Hall phase has both quasi-particle and quasi- 
hole excitations with the correct charge and statistics, and [ii) the gap towards the quasi-particle and the quasi-hole 
excitations vanishes at the critical point. 



II. THE RENORMALIZATION GROUP FORMULATION OF THE MODEL 



We now proceed with a renormalization group analysis of S. Simple power counting shows that the Chern-Sinions, 
Coulomb interactionsE3 and all three kinds of disorders are marginal at tree level in d = 2, and so loop expansions 
are required and useful. Power counting also shows that a short-range four-fermion interaction term is irrelevant and 
has therefore been iiealected in S; this makes the fermionic formulation of the anyon problem much simpler than its 
bosonic counterparttSEZlE^I. 

We assume all the three kinds of disorder satisfy Gaussian distribution with zero mean and variances Am, Ay, A^. 

< M{x)M{x') > = 1^m5'^{x - x') 
<V{x)V{x') > ^Av5'^ix-x') 

< A,{x)Aj{x') > = SijAA6'^{x ~ x') (2) 

By introducing replica a,b — 1, 2, • • • , n and doing quenched average over the above Gaussian distributions, we get 



d'^xdr 

(Pk duj 

4^2^ 






ika^{~k, -uj)a1{k, w) + -a^i~k, -wX(fc, w) 



d'^xdrdr' 
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+ Aa^<^[V'„,(x, r)7,;V'^(a:, r)][i/'„(x, T')7»V'm(a;, r')] 



(3) 



The loop expansion requires counterterms to account for ultraviolet divergences in momentum integrals; we write 
the counter terms as 



ScT —Id xdr 



a{Za - l)V'l7o9oV'm + {Z2 - l)V'l7»<9jV'm 



-^{Zf - l)g^^/^ai/^ao>:7oV'™ - ^{Zl l)g///2ai/2<V:7.C 



d'^xdrdr' 



[Zm - l)AMM'[^™(a^,T)^^(x,T)][V'„(x,r')<(a:,T')] 

- [Zv - l)Ay^'^[i/'°,(x,r)7oV'™(a;,T)][V'm(a;,'r')7oV'!'„(a;,T')] 
+ [Za - l)^A^i'WL{x, r)7, V™(x, r)] [ip^'^ix, r')7^^™(a:, r')] 



(4) 



In general, counter terms for the last two gauge field terms in S should also be considered. However, it was shownE2l 
that at least to two loops, there no divergences associated with these terms. The Ward identities following from gauge 
invariance dictate 



Zi — Za, Z^ — Z2 

Using these identities, we relate the bare fields and couplings in S to the renormalized quantities by 

i'mB = Z2 1pm 

as = {Za/Z2)a 
qB^qiJ^^'^Za/Z^Y'^ 
gB^9ii'/HZ2/ZaY'^ 
^MB — /^'^Am-^m/^2 

Ave = ii'AvZvlZl 
Aab = ^i'AAZA/Z| 



(5) 



(6) 



Notice that these relations imply that for the statistical angle 9/N = qg/N we have 9b = dfJ.'^ even in the presence 
of the Coulomb interaction and the disorders; so in d = 2 this angle is a renormalization group invariant, which is 
expected on general physical grounds. j— , 

The dynamic critical exponent, z is related to the renormalization of a byo 

d . d Z2 

z = 1 — fJ--r- In a = 1 — fi-r- In — (7) 

dfi dfi Za 

We will find it convenient to express the loop expansion in terms of the "fine structure" constant w = g^/lG, and a 
central object of study shall be its /?- function f3{w) — ^{dw / dji) . By comparing (R) with relationships between bare 
and renormalized quantities quoted above we see that 

z = l- P{w)/w. (8) 

Finally, the critical exponent v is related to the anomalous dimension of the composite operator ^ip by t^^^ — 1 = 
/^(dlnZ^ , /d/i); the renormalization constant Z-r. can be calculated by inserting the operator into the self-energy 
diagrams. 

III. ONE-LOOP CALCULATION 

We begin the explicit calculation of the renormalization constants by considering a direct perturbative expansion 
in the Coulomb fine structure constant w, the statistical angle 9 and the three kinds of disorders Am, Ay, A^. At 
one-loop order (Fig. 13, Fig.4a,b,c), we find no dependence on 9\ the values of the renormalization constants upto 
terms of order nP' , 0^, w9 and A^ are 

Z2 = 1 - 2w/NTTe 

Za = l-— (Am + A V + 2 A a) (9) 

7re 

We also explicitly verified that the Ward identities Eq.g hold. 

From Fig.5-Fig.l0 and Fig.l3a-d, we can calculate the three renormalization constants for the disorders ( the details 
are given in the next section where we perform a similar calculations in large N limit ): 

Zm = 1 + ^-r— (2A|f + 2AmAv - 4AmAa - 4.AvAa) - ^ 
tteAm iVTre 

Zy = 1 -—(Am Ay -I- 2AmAa + Ay Ay + 2AvAa) 

TreAy 

2 Aw 

Za = 1 ^AMAy - -— (10) 



From Eq.|6|, we obtain the four (3 functions: 



Piw) = /3P(u;) - -(Am + Ay + 2A^) 

TT 

/3(Am) - -2AM(i^^' - 1) + -{All + Am Ay - 2AmA^ - 2AyA^) 
/3(Ay) = 2Ay^^Aif!Z _ -(Am Ay + 2AmA^ + A^. + 2AvAa) 

W TT 

P{Aa) = --AMAy (11) 

TT 

Where the /?- function of the Coulomb coupling in the pure system pP{w) isEEl: 

27/)2 

pP{w)^^ + 0{w\w'9') (12) 

While the pure critical exponent i^p isE3: 

j/p = 1 - 2w/NTr (13) 



From Eq.O, we can identify the disordered critical exponent i^: 



.-1=.-!-^-^ + ^ (14) 



From the first and the third equations in EqJT^, we find /3(Ay) and P{w) are related by 

/3(Ay) = 2Ay^M _ IamAa (15) 

W TT 

In Eq.O, isetting Coulomb interaction and two of the random potentials to be zero, we reproduce the results of 
Ludwig et.aB. 

With the presence of Coulomb interaction w 7^ 0, we discuss the three cases separately: 

The system flows to a line of stable fixed points given by Am = ^- The flow trajectory is given by Am = -t-, C 



is an arbitrary constant (Fig. 16a). From Eqs,13, 14, we find v — 1 along this line, it saturates the bound i> > 2/d ( d 
is the special dimension ). 
Ay^O 

There exists a line of unstable fixed points given by Ay = 2h_ xhe flow trajectory is given by Ay — Cw^ (Fig. 16b). 
Below this line, the system flows to the origin, above this line, the system flows to strong coupling regime. From 
Eqs.|l3|, Q we find ly = 1 along this line. 

Aa^O 

The system flows to a line of stable fixed points given by A^i — ^. The flow trajectory is given by A^i — C 

(Fig. 16c). From Eqsjl^, |lj, we find v = 1 — ^^^ which continuously changes along this line. It should be noted that 
V <\ along this line. 

From EqlH it is easy to see that z = 1 on all these line of fixed points. It is also worth noting that, despite the 
value z = 1, the critical correlators are not Lorentz invariant. 

In the generic case, we do not find any perturbative accessible stable fixed point. The system flows to the strong 
coupling regime. In the next section, we resort to the large N expansion to explore the strong coupling regime. 

IV. LARGE N EXPANSION 

To understand larger values of 0, and to explore the consequences of a possible interference between the Coulomb 
interactions, disorders and the Chern-Simons term we found it convenient to perform a 1/A^ expansion. This is 
technically simpler than a perturbative two-loop extension of the computation of the last section and also automatically 
includes the dynamic screening of the gauge field propagator by the fcrmion polarizationE3. Alternatively stated, the 
so-called RPA approximation becomes exact at A^ = 00, and 1/N corrections require gauge field propagators which 
have the RPA form 

The inverse of the above matrix leads to the propagators of the gauge fields: 

Goiiq,!^) 



1 


9+f6 


Vq^ + y'' 


q^ 


I3 


1 



r l + (4)' + 



-16^ 16 ^q2+^2 

Goi{q,i^) = -Gio{q,v) 

To the order of (-i.)°A, the summation of Fig.l (a),(b),(c),(d),(e) leads to effective random scalar potential Ay: 

Av = ^^±^ (18) 

(A + w)^ ^ ' 



where (j) = {O/l&f and A = 1 + 0. 

Note that both Ay and A^ contribute to the effective random scalar potential with the same sign. 

To the order of (■^)°A, the summation of Fig. 2 (a),(b),(c),(d),(e) leads to effective random vector potential A^((5i 
^^i) (see the appendix): 






{\ + wy 



(19) 



Note that both Aa and Ay contribute to the effective random vector potential with the same sign. 

To the order of (■^)°A, the summation of Fig. 3 (a),(b),(c),(d) leads to effective random Scalar- Vector potential 

^ q 

^ ^^ [-Av + Aa{1 + w)] (20) 



" [X + wf 

It should be noted that both A^ and Ay contribute to the effective random SV potential with the opposite sign 
and there is no such random SV potential in the original action Eq.pl Due to the interference between disorder and 
CS interaction, this potential is generated. But it can only appear as internal lines. Putting = (namely, no CS 
interaction), this potential vanishes. 

These three effective potentials plus the random mass potential and the four RPA gauge fields propagators in Eq.[l7| 
are the building blocks in the following Feymann diagrams. Note that, however, only hare disorder potentials can 
appear as external lines of any Feymann diagram. 

In the following, we discuss N = oo case first, then we discuss 1/iV correction. It is well known that only primitive 
divergences of Feymann diagrams are needed in the RG flow equation. 

A. N = oo limit 

Fig. 4 is the contributions to the Dirac fermion self-energy from the effective random potentials. The divergent parts 
from the effective random SV potentials (d) and (c) vanish. Wc find the two constants defined in Eq.0: 

Zo, ^ 1 - —{Am + Av + 2Aa) (21) 

7re 

In the following, the notation A ^ B means that A rcnormalizcs B. 

Fig. 5 is the renormalization from the random mass ( In this figure and the following figures, we do not draw 
explicitly the diagrams with the exchange of leg3 and leg 4) . 



5a= 


7re 


-^A^ 


5c= 


M-^^« 


■^ Am 



(22) 

71 t 
^ A 2 

In all. Fig. 5 contributes — — -^ Am- 

Fig. 6 is the renormalization from the interference of random mass and effective random scalar potential. Note the 
effective random scalar potential can only appear in the internal loops. 

6a— 6b — 6c — 6d = > Aa 

ne 

2Am~Av 

6e= 6/ = > Am 

Tie 

Qg=Qh=^^^^^Av (23) 

In all, Fig.6 contributes ^AaiAv _^ ^ 4AmAv _^ ^ _^4AmAv _^ ^ 

Fig. 7 is the renormalization from the interference of random mass and effective random vector potential. 



2AmAa . 2AmAa . 

7a= 7b = > Am, > Ay 

7re ne 

2AmAa . 2AmAa . 

~c= "d = > Am, > Av 

Tie ire 

4AmAa . 

7e= 7/ = > Am 

ne 

7g= 7h = {) (24) 



In all, Fig.7 contributes ^^AmAa. ^ Am, ^^"^-^ -> Ay. 
Fig. 8 is the renormalization from the random scalar potential. 



Ba= -86 = — 5^ ^ A 



A2 



7re 



8c^8d=^^^^^Ay (25) 

7re 



In all. Fig. 8 contributes 



4AvAv 



Fig. 9 is the renormalization from the interference of effective random scalar potential and effective random vector 
potential. 

2AyA^ 2AyA^ 

9a= 9o = > Am, > Ay 

7re Tre 

Q Q^ 2AyAA . 2AyAA . 

9c= 9a = > Am, > Ay 

7re 7re 

4AyA4 

9e= 9/ = ^--^ ^ Ay 

7re 

9g= 9/i = (26) 

In all, Fig.9 contributes -^AvAa _^ ^ sAvAa _^ ^ 

Fig. 10 is the renormalization from the effective random vector potential. 

4A2 
10a= -106= ^ -^ Aa 

7re 

10c=10d=0 (27) 

In all. Fig. 10 does not contribute. 

Fig. 11 is the renormalization from the effective random SV potential. 

lla= llrf = -lie = -llh = ^ -^ Aa 

A2 

116= lie = 11/= 11.9 = -— ^^ Am (28) 

27re 

2 X 2 

In all. Fig. 11 contributes — -^ Am- It can be shown that only when SV potential get paired, there are non-zero 

contributions. 

Fig. 12 is the renormalization to the random mass from random effective scalar, vector and scalar-vector potentials. 

12a=m=i^^^^^AM 
7re X + w 

126= 12c = 12/ = 12g = ^^mAc V^ ^ ^^ 

7re X + w 
I2d= Ue = _^(^^MAy J ^ ^^ ^^g^ 

7re X + w 

In aU, Fig. 12 contributes ~^^^[^Av - {w + </))A^ - V^Ac] ^ Am- 

7 



From Fig5-Fig.l2, we find the three constants Zm, Zy . Za defined in Eq.^: 



Zm = 1 + ^-r— (2A^f + 2AAf Ay - 4AmAa - 4AyAA - A^ - 16AmAk) 
tteAm 

Zy = \ -— (AmAv + 2AmAa + AyAy + 2A,/Aa) 

TreAy 
Za = \ ^AAf Ay (30) 



TreA^ 



Where Kk is given by 



~ Ay0(3 -w-^) + ^aW-{w + (/>)(! + wf - 2^(1 + vS)\ 

(A + wy 



From Eq.|6|, we can find the following /3 functions at N — oo. 



P(^) ^ --(Am + P{Av) + 2P(A^)) 

TT 

/3(Am) = -Ai, + Ia^, + lAyA^ + -(Am - 2A^)P(Ay) - -(Am + Ay)P(A^) 

TT TT TT IT TT 

2 - ~ Ifi 

- -AcP(Ac) - -AmP(Aa') 

TT TT 

/3(Ay) = --Am Ay - -AyP(Ay) - 1(Am + Ay)P(A^) 
TT IT TT 

P{Aa) = --AMP(Ay) (32) 

TT 

Where the function P(A) = Ay^ + A^^ + «;|| + 2(/.||. 

We shall discuss the implications of this equation after considering 1/A^ corrections in the next subsection. 

B. 1/A'^ correction 

In this section, we consider 1/A^ correction to A^ = oo results. The small parameters are 1/A^, Am, Ay, A^i. We 
expect that if there are fixed points, the fixed points values of A'** are of the order 1/iV. 



Fig. 13 is the contribution to the Dirac fermion self energjj-Jrom 1/A^ fluctuation of gauge fields given by Eq, 17 
Actually, Fig. 13b and Fig. 13c are convergent. The results arec^l 

*A 1 /2u; IGw^A e^C e^E\ 

^ ~ TVTre ^ NTTe \\ tt\ IGtt IGtt ) 

"" iVTre" NTTe\ ttX Wtt '^ Wtt J ' ^^^' 

where A = 1 + (6'/16)^ and the functions A,B,C,D,E~A + B,F — B are given by the formal expressions 

f\ 4x2(1^) /-i (l-a;^)(l-6:r^+:r4) ^^ 

Jo (l + a:^)'^ Jo [l + x^y 

f^ 4x2 fi (I- 6x'^ +x'^) 
C= dx- -^--—f{x;w,e), D= dx-^^ 2^n L f(^:^^(^) (34) 







' (1 - x2)(l + a;2) ' "■ ' ' " J„ ' (1 - x2)(i + 2;2) 



with f{x; w, 9) — (A(l + x"^) + w{l — x"^))^^, and the variable x represents an intermediate frequency. Note the two 
constants C, D are divergent: this divergence is due to the singular effect of frequencies |w| :s> k. However, as shown in 
Refo, these divergences are gauge artifacts and cancel in the /3-function and in any physical gauge-invariant quantity 
like v, z or Uy . The divergences however do infect the anomalous dimension of the field operator ip: this is as expected 
as the propagator of ip is clearly gauge-dependent. 

Fig. 14 is the renormalization to Am from the 1/A^ fiuctuation of gauge fields. Actually, Fig. 14b and Fig. 14c are 
convergent, the divergent parts are 



^ 



M 



Hire 
Where the function G is given by 



J_ 4u; Aw{2A + B) 9^{2C + D) e^{2A + B) O^G 



G= / da;[((/)-l)(l-.T2) 

JO 



1 + a;^ 



(35) 



(36) 



It is easy to see Fig. 14 is exactly the same diagram for the calculation of Z^^ as expectedE2l, therefore Z2Z^^ = 
1 - *A//iV7re. 

In Figl4, replacing Am line by Ay and A^ lines, we can repeat the same calculation. In fact, the divergent parts 
should be ^^ and -^^ respectively as dictated by Ward identities. It can be shown exphcitly that Fig. 15a + Figl5b 
vanishes as dictated by Ward Identities. 

Adding the 1/N correction to the renormalization constants calculated in the last section, we obtain 



Nne 
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Zc. = l (Am 

Zm = 1 + ^^(2A 
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Za = 1- 



ttcAm 
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TreAy 
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Nne 
2AmAv - 4AmAa - AAvAa 



{AmAv + 2AmAa + AyAy + 2AvAa) 



Yc - 16AmAk) 

2£y 

Nne 



2* 



M 



Nne 



-AmA 



2^A 



M^y 



tteAa "' ' Nile 
From Eq.y, we can find 1/A^ corrections to the j3 functions in Eg. 32 



(37) 



I3{w) = /3P(u;) - -(Am + P{Av) + 2P(A^)) 

TT 

/3(Am) = -2Am{v-^ - 1) + -A2, + -Al + -AyA^ + -(Am - 2A^)F(Ay) 

TT TT TT TT 

4 __2~~1fi 

- -(Am + Ay)P(A^) - -AcP(Ac) - -AmP{Ak) 

TT TT TT 

BP(vi) 2 2-4 

/3(Ay) = 2Av^i^^ - -AMAy - -AyP(Ay) - -(Am 
W TT TT TT 

P{Aa) = --AmP{Av) 



Av)P(A^) 



(38) 



The P'P{w) in Eqj3^ is the /3-function of the Coulomb coupling in the pure case which was calculated in Ref.Ej: 



f3P{w) 



2w^{l 



Ntt^X^ 



32w(i 



TT — 16w / dx 
Jo 

1 , (l~a;2)(-l + 10a;2 
ax 



l-x' 



1 



A(l + a:^) + f(l-x^) 
(A(l+x2)+u;(l-a;2))2 



1) (l + a;2) + f(l-.T2) 



Ntt^ Jo (1 + x2)3 

The pure exponent i^p in Eq.pq is given by 

2w(l - (j)) 



1 



Ntt^X^ 



TT — 16w I dx 





(A(l+a;2)+w(l-a;2))2 



l^X^Y A(l+x2) + f(l-x2) 



l + a;2y (A(l + x2)+'u;(l-a;2)) 



2^^2 



320 /■ , .1 



l + 2;2 + f(l-2;2) 



1920 
5120(1 



dx{ 



1 — x 



"^^^(A(l + a;2) + w(l - a;2))2 
2 l + a;2 + H(i_a;2) 



l + a;2'(A(l+2;2)+u;(l-.T2))2 



dx 



(l-a;2)(l + a;2) 



7V7r2 7o (A(l + a;2) + w(l-x2))3 



(39) 



(40) 



The explicit expressions for P"^ are given by 

P(Ay) = — ^-— [Ay(l -w-i4>) + Aa<^(3 + W-^)] 

P{Ka) = — ^-^[Ay</>(3 + w-4>) + Aa(1 + w)[{\ + wf - (3 + w)4>]] 

P{Kc) - (/^)3 [-2Ay(l - (A) + A^[(A + u;)«; + 2(1 + u;)(l - ^)]] 

F(Ak) = 7-^— ^{Ay(/)(9 - 4w - 50) + Aa[5(/)2 - (1 + w)2(2u; + 3(/.) - 60(1 + u;) - 2w{w + (/>)(1 + w) - 20u.]} 

(A + wy 

- y^^J)? {Ay0(3 - ^ - 0) + A^[02 _ ei„ + 0)(1 + u;)2 - 20(1 + w)]} (41) 

From the second equation in Eq.pq, we can identify the disordered critical exponent 

^.-1 , ^,-1 Am P(Av) ^ 2P(A^) ^ 8F(Ak) ^^2) 

TT TT TT TT 

At fixed points, substituting ^ ^' — ^(Am + i^(Ay) + 2P(A^)) into Eq.En, we can simplify the above equation to 

.-^^.-^^'-^^M^^I^M (43) 

^ TT TT 

Where u'^ is listed in Eq.l2 of Ref.B 

_i 1280 /"i (l-x2)(l + 6a;2+x4) 1 + x^ + H(i _ a;2) 



!>"' - 1 = -^^ / dx 



P Nt:^ Jo (l + a;2)3 {X{1 + x^) + w{l ~ x^))^ 



5120(1-0) r, (l-^-)(l+^-) (44) 



7V7r2 7o (A(l+a:2)+w(l-a;2))3 

It can be checked that Eq.^ should be replaced by: 

d(w) 4 
/3(Ay) = 2Av^^^ - -AmP(A^) (45) 

W TT 

In Eq.p^, we expand P{w) to order 1/N and A, expand /?(A) to (l/iV)A, A^. It is easy to see there should be 
no {1/N f terms in /3(A), because interactions do not generate disorder (or equivalently disorder do not generate 
interactions ). Note that the small parameter 1/N plays a similar role to the small parameter e in the conventional e 
expansion, namely, we are trying to locate the fixed points at A'^ at the order of 1/iV. 

We now turn to the physical implications of our main results (|38| ) , (|39|) and ( [40[ ) . 

If there exists only random mass, namely. Ay — A^ — 0, Eq.|38| simplifies to 

P{w) = f3P{w) - -Am 

TT 

I3{Am) = -2AM(i^„"' - 1) + -All = -^^AiiK' - 1) (46) 

^ TT ^ 

Setting all the disorders vanishing, the authors in RefHO found a line of pure fixed points given by f3^{w) = 0. 
They also found that Vp > 1 along the fixed line, therefore concluded that this line is stable against weak random 
mass disorder from Harris criterion. Here, we explicitly write down j3{Am) to the order (l/iV)Ajv/, A^ and reach a 
stronger statement that there are no other fixed points except this line in the weak coupling regime. 

Comparing Eq.^to Eq.|38|, we find that there is no 1/N correction to j3{Aa), therefore A^ is always marginal. The 
1/N correction to /?(Ay ) is simply pP{'w)/w, therefore Ay is marginal along this line, irrelevant(relevant) above(below) 
this line. Actually, these results are expected from Ward identities. We also find that 1/N correction to /3(Am) is 
just Vp^ — 1, which is consistent with Harris criterion. 

With the presence of all the three disorders, we expect that the pure fixed line is unstable. 
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1. Integer Quantum Hall Transition ((j> = 0^ 



First, considering 
f3P{w) is simplified t' 



transition out of the integer quantum Hall state, 9 — 0, which implies = 0, A = 1. The 



pP{w) 



2w 



IQw I dx 



l+a;2 + |(l_a;2) 

(1 + a;2 + w(l - x'^)y 



(47) 



Eq.BO reduces to i'^ — 1 = f3P(w)/w. The simple analysis of (|39|) shows that pP{w) > for all w > 0; for small 



w we have [3p{w)/w = 2w/{N'k), in agreement with one-loop result ([12[), while for w ^ 1, /3(u') = A/{NTrw). So the 
only pure fixed point remains at w = 0. The whole picture of (3p(w)/w is drawn in Fig. 17a, c, it increases linearly first, 
reaches a maximum value 0.20/N at w = 1.31 and eventually decays as l/w. 

It is easy to see that Ay = jj^,P{Av) = Avjy^]Aa = P(Aa) = Aa;Ac = P(Ac) = 0;A;f = 
-AyiY^ji^, P{Ak) — —AA ^l^^,y] ■ Substituting these expressions into Eq.|3^. we find that Eq.|3^ is simplified to 



' l+io' 



f3iw) 



pP{w) 



1 



(A 



1 



P{Am) = -2Am 
p{Av) = 2A^ 



IT 

(3P{w) 
w 
pP{w) 



M 



(1 + W)3 



-A^ 

TT 



M + -(Am 

TT 



2A, 



2A4)A 



y 



/3(Aa) = — AmAv- 

TT [i 



' -AmAv 

TT 
1 — W 



W 



(1 



(A 



M 



4 4 

+ -AmA^(3- 

TT (1 + W)^ 

Av)A4 



w)^ 



Eq.Ea is simplified to 



Eq.fla is simplified to 






TT 



2A, 



TT (1 + w)'^ TT 



8Aa wiw + 2 
n (l+w)2 



/3(Ay) = 2Ay^^ - -AmAa 

W IT 



We discuss the three cases separately: 

Am ^ 

The system flows to a line of stable fixed points given by Am 
fiow trajectory is given by Am = -^, C is an arbitrary constant 
(namely, independent of large N limit ). 

This line is unstable against smallAy and A^i. 

Ay^O 

There exists a line of fixed points given by Ay — tt ^ , _ ^ ^' 



.i^M 



Like the one-loop result. 



1/iV (Fig. 17a). We suspect that ly 



(48) 



(49) 



(50) 



1 and the 
1 is exact 



which approaches oo as w 



1 . Like the one loop 
1 is exact 



result, ly — 1 and the flow trajectory is given by Ay — Cw-^ (Fig. 17b). Again, we suspect that u 

The lower part of this line (thin part) is unstable, the higher part of this line (thick part) is stable. The system 

either flows to the origin or flows to the higher part depending on the initial condition. 
This line is unstable against small Am and A^. 
A^/0 

The system flows to a line of fixed points given by Aa = f ^ . Like one loop result, the flow trajectory is given 
by Aa = C (Fig.l7c). However, unlike one loop result, v^^ = 1 - ^Aa(i _ ,^^^.i ), li w > \/2 - I ^ 0.41, v > I. 
The left part of this line (thick part) is stable, the right part of this line (thin part) is unstable. At weak disorder, 
the system cither flows to the left part of the line or to strong coupling regime depending on the initial condition. At 
strong disorder, the system always flows to strong coupling regime. 

This line is stable against small Am and Ay in the range 1 <rW < 1.31. This stable region may control the integer 
quantum Hall transitions observed in real experimental systemall. 
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As shown first by Ludwig et aJ3, the random gauge fixed fine is unstable against Am and Ay ( see also Eq.|TT| ). 
Due to the Coulomb interaction, we find there is a small part of the fixed line 1 < w < 1.31 which is stable against 
small Am and Ay. 

From Eq.^, it is easy to see that z = 1 on all these line of fixed points. 

Unfortunately, we are still unable to find a generic fixed points with all the couplings non-vanishing. These generic 
fixed points may be either unaccessible to the method developed in this paper, or simply do not exist in the real 
experimental system. 

2. Fractional Quantum Hall Transition ((f> > 0^ 

Turning to the fractional case with a non-zero 9, we start with the simplest case w — Ay = A^ = 0, /3(Am) 
simplifies to 

P{Am) = -2Am{K^ - 1) + -Al^ (51) 

TT 

Where Vp is the exponent in the absence of Coulomb interactions {w — 0): 

512</)(l-20) 

When < (f) < 1/2, Vp < 1, the pure fixed point is unstable, the system flows to a line of fixed points given by 

^^^ ^ jV37rA3 • Fi'om Eq.wand the fact ^ \w=o = 0, we obtain z — 1 + Am/tt > 1 which continuously changes 
along this line (Fig. 18). From Eq.^, we get v = 1 along this line. When > 1/2, the pure line of fixed points with 
z = 1 is stable. 

It is easy to see that this fixed line is unstable against small (w, Ay, A^). 

Unlike the Coulomb interaction case ( = ), if Ay ^ or A^ ^ 0, then all the three disorders are generated, 
this can be easily realized from Eq.^. So we have to investigate the generic fixed points of Eq.|3^. 

From Eq.pSl I3{Aa) — implies either P(Ay) = or Am = 0. In the following, we discuss the two cases separately. 

P(Ay) = 



From Eqkg and Eq.38, I3{w) = /3(Ay) = implies that P{Aa) = 0, (3'p{w)/w = Am/tt. Finally, (3{Am) = implies 



2AM{f>-' - 1) + -~Al + -AyA^ - -AcP(Ac) - -AmP{Ak) = (53) 

TT TT TT TT 



The disordered critical exponent Eqji^ is simplified to 



^ TT 

It is easy to see that P(Ay) — P(Aa) = implies the following equation: 

^ ^ Ay ^ cj,{3 + w-^) ^ {l + w)[{3 + w)<j)^{l + wf] 

Aa (11 + 30-1 4>{3 + w-4>) ^ ' 

Namely, x = 1 + w should satisfy the fourth order equation: 

x'^ + 2(0 - l)x^ - 2(j?x^ - 20(0 - 1)(0 + 2)x + 0^(2 - 0)^ = (56) 

with the constraints 3 + w > 0, > ^^^, > „ i™ . 

From Ref.El, pP{w)/w = Am/tt > imphes that < 2, therefore w < V^. 

If = 0, Eq.|5^ reduces to x^{x — 2) = 0, therefore tu = 1. If < < 1, there is no real root which satisfies both 
X > 1 and the constraints. 

If = 1, Eq.p9 reduces to {x^ — 1)^ = which implies that w — 0. Substituting (0 = 1, w = 0, t = 1) into Eq,41, 
we find Ac = -P(Ac) = P(Ak) = 0. 

If 1 < < 2, there is only one real root with x > 1. We also find P{Ax) < 0, therefore ly > Dp > 1 in this regime. 

Unfortunately, when substituting the (0, w, t) into Eq.pSl we find that the left hand side of the equation is always 
positive. Therefore, we conclude there is no perturbatively accessible fixed points with Aj\/, Ay, A^ > 0. 
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From Eq.f45|, we see that P{w) — implies /3(Ay) = 0. From Eq.^, it leads to 



dP(w) I ~ 

^ ^'^> =-(P{Av) + 2P{Aa)) (57) 



W TT 

From Eq.^, ^{Am) — leads to 

-A^, + -AvAa - -AaPCAv) - -AvPCAa) - -AcP(Ac) = (58) 

TT TT TT TT TT 

From the last equation, given (</>, i), we can determine w, then substituting {(p^w^t) to Eq.pTJ, we can determine 
A^, Ay . In the following, we consider the two end lines separately: 

(a) Ay = ( namely t = ) 

If = 0, Eq.M becomes an identity, we recover the results of the IQHT (Fig. 17c). 

If </) > 0, the solution of Eq.|| is w ^ \[icj) + ^(30)2 + 4(1 - </)) - 4] > (namely, cj) > 3/5 ). Substituting this 
expression into Eq.pTL we can determine Aa and the constraint (p < 1-3. Thus cj) must satisfy 3/5 < (f> < 1.3. From 
Eq.[43|, we find this line is stable against small Aj\/. 

(b) Aa = (namely t = oo ) 

If = 0, Eq.^ becomes an identity, we recover the results of the IQHT (Fig. 17b). 

If (/) > 0, the solution of Eq.pSl is w = 3 — 50. Substituting w = 3 ~ 5(j) into Eq.|57|, we can determine Ay: 

^-^ = i7(fz0yT- Therefore must satisfy 1/6 < < 3/5. Eq.|| becomes i^-'^ = i>-i + ^f/ (i-^)^ > ''p^ > 1> 
therefore, this line is unstable against small Am- 

We conjecture that there is a fixed plane which connect the above two end lines at t = and at i = oo (Fig. 19). 
The shaded (unshaded) region has i/ > l{i' < I) , therefore is stable (unstable) against small Am- Numerical analysis 
is needed to determine its precise boundary. The stable region may control the fractional quantum hall transitions 
observed in real experimental systemad. 

V. CONCLUSION 

Recent experiments indicated that the transitions between two Quantum Hall states or between a quantum Hall 
state and an insulating state may be described by quantum critical theories. In these theories, different FQH states 
and insulating states are considered as different ground states of the electron systems. The three important questions 
remain unsolved on the nature of these quantum phase transitions are: 

(a) The effects of the quasi-particle statistics 

(b) The effects of long-ranged Coulomb interaction on the transitions 

(c) The effects of all kinds of disorders l— , 

Answering the three questions at the same time seems a forbidding task at this moment. RefH3 investigated the 
combined effects of (a) and (b) in a Dirac fermion model and found a line of fixed points. Along this line, both 
Chern-Simon interaction and Coulomb interaction are non- vanishing, the dynamic exponent z = 1. In this paper, we 
make a serious attempt to study the combined effects of (a), (b) and (c) in the Dirac fermion model. We perform a 
renormalization Group analysis by the systematic perturbative expansions in 1/7V ( A^ is the No. of species of Dirac 
fermions) and the variances of three disorders Am, Ay, A^. we find that Am is irrelevant along this line; there is no 
1/N correction to I3{Aa), therefore A^ is always marginal; Ay is marginal along this line, , irrelevant above this line, 
relevant below this line. With the presence of all the three disorders, the pure fixed line is unstable. 

In IQHT, in the three special cases, we find the three non-trivial lines of fixed points with dynamic exponent z — 1. 

The fixed line in (Am, w) plane has v —\ and is unstable against small (Ay, Am)- 

The fixed line in (Ay, w) plane has v — \ and is unstable against small (Am, Aa)- 

Most interestingly, the fixed line in (A^, w) plane has continuously changing v and is stable against small (Am, Ay) 
in the small range 1 < w < 1-31 (Fig. 17), this stable region may control the integer quantum Hall transitions observed 
in real experimental systemad. 

The results may be relevant to the IQH to insulator transitions. It may also be important to high Tc superconductors. 
Because it was well-established that high Tc superconductors-have a d-wave order parameter and its quasiparticle 
excitations are described by 2 -f 1 dimensional Dirac fermionsEJ. 

In FQHT, setting Coulomb interaction to zero and Ay — Aa = 0, we find a line of fixed points with v — 1 and 
z> \ which continuously changes along this line (see Fig. 18). This line is unstable against small (w. Ay, A^i). 
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Most interestingly, setting Aj\/ = 0, we find a fixed plane with z = 1, the part of this plane with i/ > 1 is 
stable against small Am- This stable region may control the fractional quantum Hall transitions observed in real 
experimental system^. 

Unfortunately, we are unable to find a generic fixed points with all the couplings non-vanishing. These generic fixed 
points may be either unaccessible to the method developed in this paper, or simply do not exist in the real experimental 
system. However, by looking at carefully the divergent structures of all the relevant Feymann diagrams, we show the 
model is renormalizable to the order (1/A^)A, A^, (1/iV)^ discussed in this paper; we do bring out the systematic and 
elegant structure which describes the interferences between Chern-Simon interaction. Coulomb interaction and the 
three kinds of disorders. We believe that the structure is interesting in its own right and may inspire future work to 
study this difficult problem. 
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APPENDIX A: THE PROOF OF EQUIVALENCE OF RANDOM VECTOR GAUGE POTENTIAL AaSij 

AND Aa{Sij - ^0J-). 

We can decompose the random gauge field Ai{x) in Eq.|l| into transverse and longitudinal components: 

A, ^ Af + Af = e^.d.x' + d^x"- (Al) 
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Where x^{x):X^{^) satisfy 

< X^{^)x^ix') >=< x''{x)x''{^') >= -A^ln \x - x'\ (A2) 

From the above equation, it can be shown easily 

Kikj 



< Af{k)Aj{k') > = Aa{6., -i^), 



< Afik)Afik') > =Aa'^ (A3) 

Adding the two equations above leads to 

< A,{x)A,{x') >-< Af{x)Aj{x') > + < Af{x)Afix') >= A^5,y5'^(x - x') (A4) 

which is the third equation of EqJ2[ 

By gauge transformation, Af{x) = diX^{x) can be removed, so < A^{k)A^{k') >= A^-p^ should not make any 
contribution to gauge-invariant fmantities. This fact is similar to the "running gauge fixing parameter " in usual 
relativistic quantum field theorjcj. This point can also be demonstrated in the following specific example. 

Let us evaluate its contribution to fcrmion self energy Fig. 4c, the divergent part is 

- — (7ow + 7i/ci) (A5) 

It is evident that although this contribution infect the anomalous dimension of the field operator "0, it does not 
affect the dynamic exponet z which is a gauge invariant quantity. 



Figure Captions 

FIG. 1. The effective random scalar potential Kv to the order (l/A'^)''Av and (1/A'^)°Aa. The thick dashed line is the 
effective random scalar potential. The thin dashed lines are the bare random scalar and vector potentials. The wave line is the 
gauge fields propagators. 

FIG. 2. The effective random vector potential Ka to the order (1/A'')''Aa and (l/A'')°Ay. The thick dashed line is the 
effective random vector potential. The thin dashed lines are the bare random scalar and vector potentials. 

FIG. 3. The effective random Scalar- Vector potential Ac to the order {l/N^' IS.A and (l/A'^)"Av. The thick dashed line is 
the effective random SV potential. There is no bare random SV potential. 

FIG. 4. The contribution to the Dirac fermion self-energy from random potentials. 

FIG. 5. The renormalization from the random mass. 

FIG. 6. The renormalization from the interference of random mass and effective random scalar potential. Note the effective 
random scalar potential can only appear in the internal loops 

FIG. 7. The renormalization from the interference of random mass and effective random vector potential 

FIG. 8. The renormalization from the random scalar potential 

FIG. 9. The renormalization from the interference of effective random scalar potential and effective random vector potential 
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FIG. 10. The renormalization from the effective random vector potential 

FIG. 11. The renormalization from the effective random SV potential 

FIG. 12. The Renormalization to random mass from Av, Aa, Ac. In the text, it is called Ax. Compare this figure to 
Fig.l3. 

FIG. 13. The renormalization to the Dirac fermion self-energy from gauge fields fluctuations to order 1/N 

FIG. 14. The renormalization to the random mass from gauge fields fluctuation to the order 1/A^ 

FIG. 15. The fermion bubbles which contribute to the renormalization of random scalar and vector potentials. In the figure, 
M,!^, A = 0,1, 2. 

FIG. 16. The One-Loop Renormalization Group flow of (a) random mass and Coulomb interaction (b) random scalar 
potential and Coulomb interaction (c) random gauge potential and Coulomb interaction 

FIG. 17. The Renormalization Group flow to order 1/N of (a) random mass and Coulomb interaction (b) random scalar 
potential and Coulomb interaction (c) random gauge potential and Coulomb interaction. The thick (thin) lines in (b) and (c) 
are stable{unstable ) line of fixed points. 

FIG. 18. The Renormalization Group fiow to order 1/N of random mass and Chern-Simon interaction 

FIG. 19. If Am = 0, there is a fixed plane in {(j),t — Av/Aa) plane, w and Aa and Ay are uniquely determined by (0,^), 
therefore are not shown in the figure. The shaded regime with 2 = 1, j^ > 1 is stable against small Am. 
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